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. The morphology of clusters of galaxies may be described with a set of parameters which 

"iJ ' contain information about the formation and evolutionary history of these systems. In this 

paper we present a preliminary study of the morphological parameters of a sample of 28 com- 
pact Abell clusters extracted from DPOSS data, measured with a procedure based on the use 
Q-f of the CIAO-Sherpa software, developed at the Center for Astrophysics (CfA) for X-ray data 

q analysis. The morphology of galaxy clusters is parameterized by their apparent ellipticity, po- 

S_h ■ sition angle of the major axis, centre coordinates, core radius and /3-model power law index. 

^ \ Our procedure provides estimates of these parameters (and of the related uncertainties) by 

C3 ■ simultaneously fitting them all, overcoming some of the difficulties induced by sparse data 

and low number statistics typical of this kind of analysis. The cluster parameters were fitted 
in a 6 x 6 Mpc 2 region, measuring the background in a 4 Mpc < R < 5 Mpc annulus. We 
also explore the correlations between shape and profile parameters and other cluster proper- 
ties. Our results can be summarized as follows: one third of this compact cluster sample has 
core radii smaller than 100 kpc, i.e. near the limit that our data allow us to resolve, possibly 
consistent with cusped models. The remaining clusters span a broad range of core radii up to 
~ 1500 kpc, including some apparently regular clusters with well resolved core radii. More 
than 80 per cent of this sample has ellipticity higher than 0.2. The alignment between the 
cluster and the major axis of the dominant galaxy is confirmed at a high significance level, 
while no correlation is observed with other bright cluster members. No significant correlation 
is found between cluster richness and ellipticity. Instead, cluster richness is found to corre- 
late, albeit with large scatter, with the cluster core radius. Finally, in contrast with claims in 
previous works, a flat universe seems to be favoured, and in any case is not excluded, by the 
power-law index f3 of our number density profiles. 

Key words: Galaxies: clusters: general - Galaxies: clusters: individual - observational cos- 
mology 
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1 INTRODUCTION 

In CDM cosmological models, the formation of cosmic structures 
is associated with that of dark matter halos, which form hierarchi- 
cally from the gravitational collapse of primordial density fluctu- 
ations. Unlike galaxies, clusters are dark matter dominated sys- 
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terns in which the baryons are not radially segregated. Also, un- 
like galaxies, they form by dissipationless collapse. Thus, clusters 
in principle provide a more direct probe of the primordial density 
fluctuations. 

The most accurate methods to derive the mass profile of 
galaxy clusters (and also their internal structure) require either 
strong/weak lensing analysis or, under the assumption of thermody- 
namical equilibrium of the system, the velocity dispersion or the X- 
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ray gas density and temperature profiles. Apart from the validity of 
the equilibrium assumption, the acquisition of the necessary data is 
often an expensive task for large, statistically significant samples of 
clusters. Alternatively, under the assumption that galaxies trace the 
underlying mass distribution, cluster-sized dark halos density pro- 
files can be derived from galaxy number counts. While this method 
is not exempt from systematic uncertainties, it relies on few as- 
sumptions and, unlike the other methods, can be easily applied to 
large samples o f clusters extracted fr om wide-field surveys. 

iNavarro. Frenk & White! \ 1996ft showed via extensive numer- 
ical simulations that in a CDM universe, halos of different mass 
(from dwarf galaxies to rich clusters) follow a universal cusped pro- 
file (hereafter NFW): 



(1) 



Pcrit (r/r s )(l + r/r s ) 2 ' 

The NFW model and its universality have been challenged 
by several authors, especially because of the rotation curves of 
low surface brightness galaxies and dwarf galaxies, which show 
a relatively flat density distribution (see e.g. , [MooreeUalJU 999t 
ISwaters. Madore & Trewhelld 1200(1 Ide Blok & Bosmdl2002h .As 
applied to galaxy clusters, even if the NFW mod el is often ac- 
cepte d as a good fit to mass density profiles (e.g. [Carlbe rg et alj 
I l997t iMarkevitch et al Il999t iGeller. Diaferio & Kurtzl Il999t 
iLewisiBuTte^S^tockeT 2003ll controversy still exists and has 
been recently revived by several authors (see for instance 
| Miralda-Escuddll995tlshapiro & Ilie-J2000t lwu & Chiuehll200ll: 
Isand. Treu & EllisfcOOllBiviano & Girardil2003h . 

From the number density and luminosity profiles of rich 
galaxy clusters in the ENACS (ESO Nearby Abell Cluster Survey) 
catalog, Adami et al. (1998, 2001, hereafter AMKB98, AMUS01) 
found that, if the faint cluster members are included, a core model 
of the form: 
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reproduces the galaxy number density profiles more accurately than 
cusped (NFW-like) models; the trend turned out to be even more 
evident for luminosity profiles. This result is in agreement with 
models predicting different evolutionary scenarios for the luminous 
and the faint populations of galaxies in clusters. In fact, AMUS01 
suggested that galaxy number density and luminosity profiles could 
be original ly cusped, as exp ected from simulations of dark mat- 
ter halos iNavarro. Frenk & White 1996l l 19971) and that, afterward, 
the cusp in the number density profile could be erased by merging 
events. This would produce a cusp in the bright galaxy luminosity 
profile, as well as the core observed in the faint galaxy luminosity 
profile through tidal disruption of small galaxies near the cluster 
centre. 

The cluster morphological parameters thus bear relevant in- 
formation on the cluster dynamical status and formation scenarios. 
For instance, the core radius r core is a characteristic scale-length of 
the galaxy distribution inside the cluster, and it has been shown to 
be tightly correlated with the virial radius (see lGirardi et all 19951 
hereafter G95). Other relevant information may be obtained from 
the cluster shapes once they are approximated as spheroids. Rich 
clusters have been found to show ellipticities up to 0.8, with a 
mean value of 0.4, and a correlation between ellipticity and rich- 
ness (the richer the cluster, the more spherical the shape - e.g . 
IStruble & Ftacladll994 Ide Theiie. Katgert & van Kamperjfl995l) . 
Moreover, the cluster shape distribution is directly connected to 
the cosmological parameters: a low density universe produces more 



concentrated, spherically symmetric clusters than a = 1 scenario 
(e.g. lEvrard et alll993^ l. 

In this work we present a preliminary study of the morpho- 
logical properties of a sample of 28 nearby compact Abell clus- 
ters extracted fro m the Digitized Palomar Sky Survey (DPOSS) 
jDiorgovski etalll999h . The advantages of using DPOSS data are 
mainly in the wide sky coverage (which allows an extended region 
around each cluster to be studied), the homogeneous photometric 
system of the whole data set, and g ood control of the cata logue 
comp leteness and selection criteria l Weir et alll995tlPaolillo et ail 

We point out that we do not attempt to test which model 
best describes the actual cluster profile. Instead, we adopt the 
beta model as a statistically adequate, analytically convenient de- 
scription of the galaxy distribution and investigate the correla- 
tions between shape, profile parameters and other cluster prop- 
erties. The comparison between different cluster profile models 
requires careful stacking of individu al clusters to obtain a suffi- 
cient S/N ratio (see for in stance ICarlberg. Yee & Ellingsonlll997t 
Ide Theiie & Katgeri l999). and is beyond the purpose of this work. 

The paper is structured as follows: in section §|2|we describe 
the data and the cluster sample, in section § [3] we describe the 
method used to model cluster radial profiles and to derive mor- 
phological parameters, while in sections § |4] and § [5] we present 
our results and summarize our conclusions. Throughout this paper 
we use Ho = 50, Qm = 1, = 0. We do not adopt the usual 
concordance cosmology in order to simplify the comparison with 
Paolillo et al. (2001), on which this work is partially based. Cosmo- 
logical parameters only affect our results via the conversion from 
angular to linear distances. Using a Ho = 70, Qa = 0.7, SIm = 0.3 
cosmology, distances would be smaller by approximately 20-25 % 
in our redshift range. 



2 DATA AND CLUSTER SAMPLE SELECTION 

This work is based on catalogues produced from DPOSS data 
toiorgovski et alj Il999h with the SKICAT package JWeir et alJ 
1995). We started with the initial sa mple of 80 Abell clusters de- 
scribed in Piranomonte et al. (2001), which had available calibra- 
tion frames and at least one reliable spectroscopic redshift (see 
IPaolillo et all 1200 ll) at the beginning of this work. We note that 
a much larger cluster sample extracted from DPOSS catalogs is 
curren tly availa ble and its analysis is postponed to a future work 
iGal et alliooll) . 

The photometric completeness limit was estimated for each 
cluster and in each band independently; typical values are m r ~ 
20 and m ~ 20 .5 in the Gunn-Thuan photometric system 
IPaolillo et ai]l200lt Paolillo et al., in preparation); K-corrections 
were not taken into account since, due to the small redshift range 
covered by the clusters in our sample (z < 0.27, with only 2 clus- 
ters at z > 0.2), they are negligible; furthermore, the lack of mor- 
phological information for the individual galaxies would introduce 
unnecessary ambiguities. 

For each cluster we extracted from the SKICAT catalogues all 
objects brighter than m 9 = 20 contained within a square region of 
10 x 10 Mpc 2 centered on the cluster centre (as listed in the Abell 
catalogue), detected in both the g and the r bands (the i-band was 
not used since it is much shallower than the others) and classified as 
galaxies in both filters. The large region around each cluster (10 x 
10 Mpc 2 ) was selected in order to sample the cluster galaxy distri- 
bution out to a significant distance from the cluster centre (1 Abell 
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radius = 3 Mpc with our cosmology), and to have a control field 
around each cluster wide enough to achieve a proper background 
determination. As discussed later (§ 13. H . the proper cluster fitting 
was performed on a smaller region (6x6 Mpc). 
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Knebe & Muller 2000). In particular, cluster substructure may 



generally be divided into four categories: dynamical subclumps 
residing in generally relaxed systems; young clusters in an early 
merging state; dynamically bound units in the cluster outskirts; 
and chance projections of groups. Discriminating between these 
four different cases is extremely difficult in the absence of velocity 
information. Moreover, subclustering at large distances from the 
cluster centre (> 1 Mpc) may oft en be due to fluctuations in the 
extended host supercluster (West & Bothun 1990). 

Obviously, the presence of substructures would make the mod- 
eling of the radial profile almost meaningless, since there would be 
no radial symmetry nor a well defined core radius associated with a 
central relaxed region of the system. For the purpose of this work, 
we excluded from our sample all clusters with obvious large sub- 
structures by visually inspecting the galaxy surface density maps 
(smoothed with a Gaussian filter with a = 250 kpc) and rejecting 
those with multiple isolated peaks within 3 Mpc of the cluster cen- 
ter (see figure^. 

While this is an empirical criterion, a more objective method 
is not straightforward in our case, due to the low number of galax- 
ies and the lack of spectroscopic information. Our criterion allows 
the selection of clusters that are broadly regular and single peaked 
on large scale, and therefore our sample includes several examples 
of regular clusters like A1914, A1835, A566 liuote & Tsall996l 
Maier owicz. Neumann & ReiprichlEool IZhou et alJl2003h . How- 
ever, we also have some systems which are not truly relaxed, gener- 
ally known from X-ray or dynamical analysis such as the disturbed 
morphology of A2061 <Pos tman 1 Gel ler & Huchr3l 19881) . the late 
merg er A2142 ( Buote & Tsai 1996), and the bimodal c luster 
A98 (Beer s. Geller & Huchral |l982; Krempec- Krygier & Krygien 
1995). We prefer to use a homogeneous criterion for the whole sam- 
ple without discarding individual clusters for which additional data 
are available, even if this means that we have to accept some degree 
of 'imperfection', e.g. low density tails or some level of substruc- 
ture in the cluster core. 

It may be argued that a smaller smoothing scale may reveal 
smaller scale substructures. In figure |2| we show the comparison 
of density maps smoothed with 250 kpc and 150 kpc smoothing 
scales for two representative cases. We find that, using a 150 kpc 
scale, only five clusters (A175, A655, A910, A 1661, A2177) reveal 
significant substructure that may in principle question their inclu- 
sion in the sample. Our data do not allow to study smaller scale 
structures in the cluster core since the distance between galaxies is 
on average greater than 50 kpc (i.e. on scales smaller than 150 kpc 
shot noise dominates the galaxy distribution). Note however that a 
smaller smoothing window only affects the sample selection and 
not the actual fitting (see section § 13. 1> . 

The final sample of clusters which will be the subject of the 
following analysis consists of 28 objects listed in tableQ 



3 THE DETERMINATION OF RADIAL PROFILES 
3.1 Model and statistics 

Modeling the 2D distribution of galaxies in clusters is a problematic 
task due to the small number statistics which affects these sparse 
data. The classical approach which obtains the number count ra- 
dial profile by binning and counting galaxies in concentric circular 
shells, results in the loss of relevant morphological information. 
Moreover, because of contamination from background/foreground 
objects, there is a significant decrease in the S/N ratio of the result- 
ing profiles. 

Maximum likelihood methods as initially proposed by Sarazin 
(1980), applied directly to the galaxy positions without any bin- 
ning, allow these problems to be partially overcome. In general, 
these methods rely on three fundamental assumptions: i) the ob- 
served positions of galaxies are statistically independent (i.e. no 
substructure is present); ii) the local background is assumed to be 
uniform; iii) the shape of the galaxy number density function is 
assumed to be known and characterized by a small number of pa- 
rameters. 

In the literature, radial profiles of galaxy clusters have mainly 
been determined by means of different implementations of this 
method. However, even if in principle it would be possible to fit 
all the parameters simultaneously, in practice the minimization of- 
ten has to be accomplished insteps, fitting the shape (centre coor- 
dinates, ellipticity and position angle) and the profile (core radius, 
power law index etc.) p arameters separately (e.g . AMKB98; see 
also G95, AMUS01 and lGirardi & Mezzettill200ll) . Furthermore, 
in such cases the strong correlation between the model parameters 
may produce spurious results. 

Taking advantage of the formal analogies existing between X- 
ray data (photon hit positions and photon energies) and galaxy cat- 
alogues (galaxy positions and magnitudes) in terms of sparse dis- 
tribution and low statistics, we tailored to our specific needs the 
Sherpa package present in the CIAO 1 software, developed at the 
Chandra X-ray Center. Our procedure allows us to fit all the pa- 
rameters of the galaxy distribution at the same time, instead of re- 
quiring a multi-step approach. 

The Sherpa software works on binned data. To minimize the 
loss of information, we binned the galaxy positions using a square 
grid of bin size 50 kpc, in order to have at most one galaxy per bin, 
except for a few clusters in which the highest density bin may con- 
tain two galaxies. This is in principle equivalent to using unbinned 
data. 

We assume that galaxy clusters can be described by a standard 
two dimensional /3-model of the form: 



fir) = f(x,y) 



E 



(1 + (r/r core )Y 



where: 



~{x,y) = 



^Jx 2 {l -e) 2 + y 2 
1 - e 

x — (x — xo) cos 8 + (y — yo) sin 6 
y = (y — yo) cos 8 — (x — xo) sin 6. 



(3) 



(4) 

(5) 
(6) 



This implies that the characterization of the galaxy profile 
of each cluster requires the simultaneous fit of the following 
parameters: 

1 We used CIAO version 2.2.1. Full documentation can be found at 
http://cxc.harvard.edu/ciao/ 
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i) shape parameters: centre coordinates {xo,y$), ellipticity e, and 
position angle 9; 

ii) profile parameters: core radius r CO re, power law index P, central 
density So at xo, yo), and local background density £(,fc 9 (assumed 
to be uniform ) 2 . 

As already explained in the introduction, our goal is not to 
determine which function describes best the profile of galaxy clus- 
ters. We are heavily limited by the small number of galaxies and by 
the sparse nature of the galaxy distribution, which makes it difficult 
to discriminate between different profiles for most clusters in our 
sample. As an example, in figure [3] we compare the best-fit NFW 
and /3-model profiles, derived through an azimuthally averaged uni- 
dimensional fit on binned data, for A286. Even with this simplistic 
procedure, which minimizes the number of fitted parameters, the 
two models are almost indistinguishable. In §|4]we will discuss fur- 
ther the possibility that our clusters follow a cusped distribution. 

The use of parametric methods has also been criticized, in par- 
ticular when not dealing with truly sparse data or theoretically well- 
established functional forms. This criticism arises mainly from the 
fact that projection effects can transform two intrinsically different 
spatial profiles to very similar projected profiles, and by the appli- 
cation of parametric methods to a typical case of an ill- conditioned 
probl em such as the radial profile estimation l M erritt & Tremblavl 
1994). Although in our case we are actually dealing with sparse 
data, so that a non-parametric approach would be quite difficult, 
we are aware that we are simply using a statistically convenient 
empirical formula in a consistent manner. 

In performing our fit we used a maximum likelihood ap- 
proach, where the likelihood function is the product of the individ- 
ual probabilities Pi computed for each bin i assuming that galaxy 
counts are sampled from a Poisson distribution: 
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-exp(-Pi) 



(7) 



In equation Pi is the sum of cluster and background model 
amplitudes, and TV; is the number of observed galaxy counts in bin 
i. The Cash statistic lCashll979l) : 



dnPi 



(8) 



is derived from the likelihood function by taking the logarithm of 
— £ and dropping the factorial term, which is constant in fits to the 
same dataset. Unlike the more traditional A^ 2 , the Cash statistic 
may be used regardless of the numb er of cou nts in each bin (for 
details, see the Sherca Reference Man ual200ll) . We note that, with 
poorly sample data, the Cash statistic works better than a \ 2 statis- 
tic, even if the latter is tailored to work with low number counts. 



3.2 Fitting procedure and validation 

The application of the CIAO-Sherpa software for the determination 
of cluster radial profiles was extensively tested on real and simu- 
lated data. A first problem to solve was how to estimate the galaxy 
background. Tests performed on the 10 x 10 Mpc 2 region (i.e. 



2 The external constraint due to the normalization of the model to the ob- 
served total number of galaxies, reduces the degrees of freedom of the prob- 
lem. Therefore, only three of the four profile parameters are actually inde- 
pendent. 



R < 5 Mpc) showed that, for the majority of clusters, the fit was 
able to recover within 10% the estimate of the background level as 
obtained manually, by measuring the average galaxy counts in an 
external annulus (4 Mpc < R < 5 Mpc) around the cluster. Both 
the manual and fitted estimate are reported in table|3| 

For low S/N objects, while still retrieving the correct back- 
ground, problems were encountered in fitting the other cluster pa- 
rameters in the whole 10 x 10 Mpc 2 region. This is a natural conse- 
quence of the fact that with low S/N objects, the background fluctu- 
ations undermine the cluster detection itself. Therefore, in order to 
be consistent, we decided to fix for all the clusters the background 
density to the value measured in the 4 < R < 5 Mpc annulus, and 
to fit the cluster in a smaller region of 6 x 6 Mpc 2 . We tested the 
effect of a background under- or over-estimation by 15% (~ 2a 
upper limit on the background fluctuations measured in our cata- 
logs over a 6 x 6 Mpc 2 window) of its measured value on the fitted 
parameters r core , /3, ellipticity and position angle: as expected, the 
ellipticity and position angle are scarcely affected, while there is 
a systematic effect on f3 and, to a lesser extent, on the core radius 
(see figure[4}- 

The 6x6 Mpc 2 was adopted as an optimal compromise be- 
tween the requirement to sample a relevant part of the cluster, 
possibly out to the radius where the galaxy density approaches 
the background level, and the attempt to maximize the S/N ra- 
tio of our galaxy counts. Note that this region corresponds to a 
radius R — 3 Mpc, i.e. one Abell radius in our adopted cos- 
molog y, and app roximately to the median cluster virial radius 
ICarlberg. Yee & Ellingson 1997). Furthermore this choice helps 
to prevent strong background fluctuations or nearby (3 < R < 5 
Mpc) overdensities from affecting the cluster parameter evalua- 
tion. The 4 Mpc < R < 5 Mpc annulus represents a region dis- 
tant enough not to be contaminated by the cluster itself while suffi- 
ciently nearby to be representative of the local cluster background. 
Since the use of a 6 x 6 Mpc region may include small local sub- 
structures (either due to the cluster or the background) which are 
not discarded by our sample selection criteria, we performed an ad- 
ditional fit on a smaller 3x3 Mpc 2 area. The results, which are 
discussed in section §|4] show that the fits performed on the smaller 
region, while severely affected by the lower statistics, are basically 
consistent with those performed over the 6x6 Mpc area. We show 
the dependence of the retrieved parameters on the size of the fitted 
region in figure [5] The estimates in the three different regions ap- 
pear consistent within the errors in the vast majority of cases. The 
f3 parameter shows indeed a large scatter, as this is the less sta- 
ble of the fitted parameters, due to low S/N and incomplete profile 
sampling in small apertures. 

The fit results do not depend on the bin size, provided that the 
bin is small enough to avoid any loss of information (i.e. ~ 1 galaxy 
bin - ). Furthermore, the shape parameters (centre coordinates, e 
and 9) derived from the fit are consistent within the errors with 
respect to profile (r core , f3) parameter variations (see figure|5}- 

We don't find any significant correlation of the fit result on the 
starting values: using different values within the permitted range 
(i.e. [0:3000] kpc for w e , [0:4] for/3, [0:1] fore, and [0:2tt] ford), 
the difference in the retrieved output is 5-20% for r core , 1-15% for 
0, less than 2% for e, and less than 0. 1 % for 9 (the larger differences 
quoted for r core and occourring for the clusters with very large 
core radius). 

Finally, we tested our procedure on a sample of 230 mock 
clusters. Mock clusters are generated following a /3-model profile, 
and all clusters are produced with centre coordinates in the centre of 
the mock catalog region, which in pixel units is (xo, yo) = (60,60), 
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ellipticity = 0.4 and /3 = 0.8, representative of typical values in the 
real sample, and position angle = 1. The core radius, central density 
and background density are randomly chosen in such a way that the 
total number of objects, the background density ai,k g , the cluster 
core radius, and the S/N reflect those of the real clusters. The fit on 
the mock catalogs is performed exactly in the same way as for real 
clusters: background density is fixed and all the other parameters 
are left free to float. 

The mean ellipticity retrieved is 0.33 ± 0.16, the mean j3 
is 0.83 ± 0.26, the mean centre coordinates are (60.2,60.1) ± 
(1.2, 1.5), the mean position angle is 1.3 ± 0.7. The distributions 
of the retrieved ellipticity, /3, centre coordinates and position an- 
gle are shown in figure |6| and are to be compared with the above 
reported true-values of (xo,yo) = (60,60), ellipticity = 0.4, posi- 
tion angle =1., and f3 = 0.8. In figureQwe show the retrieved vs. 
true core radius. These plots show that our fitting procedure yields 
reliable results without significant systematic effects, except possi- 
bly for the ellipticity which is slightly skewed toward lower values. 
This is a result of the background contribution which tends to az- 
imufhally smooth the galaxy distribution. This effect however does 
not affect our results, as discussed in §|4| The large dispersion in the 
P.A. distribution is mainly due to low S/N. In fact, additional sim- 
ulations with very high S/N mock clusters show that for e > 0.3 
the position angle is correctly retrieved within 10 degrees at most. 
Note that while the core radius exhibits a large relative dispersion 
at small radii, due to the difficulty of obtaining a reliable estimate 
of r core for small core radii (due to the low S/N within r core ), the 
retrieved value is in most cases representative of the true value, and 
the scatter found for mock clusters is consistent with the average 
measurement errors that we derive for the real sample. 

3.3 Fitting results 

Summarizing the discussion in section § 13.21 we performed the fit 
for all clusters in a 6 x 6 Mpc 2 square region, with a bin size of 50 
kpc and with the background density fixed to the value measured 
for each cluster in an external annulus (4 Mpc < R < 5 Mpc). 
Note that the use of the same physical area (6x6 Mpc 2 ) for all 
clusters implies that we are considering a different density thresh- 
old for each cluster. However, this does not affect our results since 
we sample the cluster out to a radius where the S/N is negligible, 
while the use of a fixed density threshold would introduce other 
effects, such as a dependence on the background fluctuations. 

The resulting parameters are listed in table [2] Error estimates 
are calculated by leaving only the highly correlated parameters free 
to vary (i.e. e — 9 and r core — fi — £o), and fixing all the others. 
Even if these are not the formally correct errors, they provide a 
meaningful uncertainty on the parameter estimates. The formal er- 
rors (where all the fit parameters are left free to vary) would often 
be unconstrained, given the low counts and the large number of 
parameters being fit. 

As a template case, we show in figure|8|the results of the fitting 
procedure for the cluster A98. In the upper left panel, the solid line 
traces the best-fit profile and the dots show the measured density 
profile. The data and the model are binned using the same set of 
annuli centered on the cluster centre defined by the fit. We point 
out that this binning is only used for visualization purposes and it is 
not involved in the determination of the best-fit parameters. In the 
upper right panel, we plot the cluster isodensity contours against 
the model. 

In the lower left and middle panels, we plot the confidence 
contours at la and 2a for the highly correlated parameters, namely 



e — 9 and r core — /3 (with Eo free to vary). Finally, in the lower 
right panel we plot the distribution of | cos(9 ga iaxies - 0duster)\ 
(solid line) and the expected histogram for a random distribution of 
the galaxy position angles (dashed line). For the other clusters in 
the sample we show only the three most significant plots, namely 
the cluster isodensity contours against the model, the fitted profile, 
and the | cos(9 ga iaxies — 0duster) \ distribution (fig.|9j. 



4 DISCUSSION 

In figure fTol we present the distribution of ellipticities (left panel), 
slopes /3 (middle panel), and core radii (right panel) derived from 
our fits. Even though the sample was selected to include clus- 
ters with an overall regular morphology, we find a wide range 
of ellipticities, with mean value e = 0.47 and dispersion a t — 
0.2; note that, given the systematic error measured by our sim- 
ulations (S 13.2k the real ellipticity distribution is expected to be 
slightly skewed toward larger elongations by Ae ~ 0.07. While 
thi s result is in fair agreement with the e distribution measured 
bv lde Theiie. Kateert & van Kamper] 1 19951) . it is in contrast with 
AMKB98, who found on average low ellipticities, consistent with 
e ~ 0. However, AMKB98 and AMUS01 were mostly interested 
in the central cluster region, in order to compare cusped and core 
profiles; as these authors point out, their study is limited to the very 
central region of the cluster (~ 500 kpc) specifically chosen to 
avoid substructures. This implies that, apart from possible aperture 
biases which could affect the result, the cluster region that they fit 
has intrinsically very small ellipticity. 

In fact, in our sample six out of the eight less regular clusters 
mentioned in § 2 have ellipticity above the mean. 

As noted in section § 13.21 in order to better evaluate possi- 
ble effects introduced by overlooked background substructures in 
the cluster outskirts, we performed an additional fit over a smaller 
3x3 Mpc 2 region. While the parameters derived in this way are 
more uncertain because of the lower statistics, and possibly biased 
by the exclusion of the cluster external region, the ellipticity distri- 
bution and the presence of the peak at e ~ 0.5 are confirmed also 
using the parameters derived from this smaller region (the dashed 
line in figure fTol . Furthermore, in figure fTTI we plot the position an- 
gles 9 derived within the 6x6 Mpc 2 region against those derived 
within the 3 x 3 Mpc 2 region (clusters A2061, A2083, A2178 and 
A910 are not included because no reliable parameters could be de- 
termined for them in the 3x3 Mpc 2 region). The two estimates are 
fully consistent and the few points scattered above and below the 
line are clusters with e ^ 0.25, i.e. objects for which the error on 
the position angle is systematically larger (the higher the ellipticity, 
the better constrained the position angle of the major axis). This 
result shows that there is no significant substructure in the clus- 
ter outskirts which would affect the cluster P.A., even though the 
contribution of cluster galaxies is non-negligible between 1.5 and 
3 Mpc from the cluster centre: a significant fraction (between 20 
and 60 per cent, with a median value of approximately 50 per cent) 
of cluster members is found in this external region. Therefore, any 
study of cluster structure sampling less than 1.5 h _1 Mpc in radius 
is missing a significant cluster contribution. 

The power-law slope /3 has a peak at ~ 0.7, with /3 = 0.8 and 
dispersion ap = 0.27; both the average value and the shape of the 
dis tribution are in very good agreement with the results obtained 
bv lPopesso et al. I 120041) based on the ROSAT-SDSS galax y clus- 
ter sample, as well as with the results of lGirardi & Mezzettil 1200 ll) 
based on a sample of moderate/high redshift clusters . On the other 
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hand, our results are again in contrast with those of AMKB98, who 
found P — 1.02±0. 08 for their sample of 60 clusters. Even though 
our distribution is quite broad, the average value of p is inconsistent 
at the > 3a level 3 and none of the clusters in the AMKB98 sample 
has P < 0.84. While our catalogues are somewhat deeper (~ 0.5 
mag) than the ones used by AMKB98, the main differences be- 
tween us and AMKB98 are that the latter sample a smaller region, 
thus excluding the contribution of the cluster outskirts, and by the 
different approach used to estimate the background. Note however, 
that the fits performed on our smaller 3x3 Mpc region, which 
is closer to the AMKB98 area (r > 5r core ~ 1000 kpc 4 ), still 
yield P = 0.75. The two P estimates for the individual clusters are 
consistent within the errors and there is no evident correlation of P 
with the size of the sampled region. Note that, according to figure 
10 of AMKB98 (based on simulations by ICrone & Evrardll994t 
Jing et aljll995l: iNavarro. Frenk & WMtefll995t IWalter & Klypir] 
1996) our P distribution is compatible with flat universe models ei- 
ther with zero (standard CDM) or non-zero cosmological constant 
(fi ro = 0.2 0.3 and A = 0.7 ~ 0.8), in contrast to the AMKB98 
result which favoured an open cosmology. We remind, however, 
that the validity of cosmological constrains obtained through this 
method is still debated, and depends on the used simulations. 

The median core radius is r core — 310 kpc, but the dis- 
tribution is strongly skewed toward low values with a marked 
peak at r core ~ 100 kpc. Approximately 1/3 of our sample has 
r C ore < 100 kpc, corresponding to 2 pixels in our fitted maps; for 
such clusters r core is thus only marginally resolved, as discussed 
in § 3.1, so that we cannot exclude a cusped profile. On the other 
hand 7 clusters have core radii greater than 500 kpc; among these 
we find six of the less regular systems discussed in § 2. 

Some correlation between the dynamical state of the cluster 
and its core radius is clearly expected, and in fact the most regular 
clusters (i.e. A79, A763, A971, A1677, A2065, A2083, A2223) all 
have very small core radius. The intermediate values of core radii 
refer to clusters which are generally single peaked, but exhibit some 
broader asymmetric overdensity region, or a filamentary structure. 
These systems also include the examples of known regular clusters 
mentioned in § 2. 

Our core radii, including the whole sample, are on average 
larger than those measured by AMKB98 and AMUS01. How- 
ever, our mean core radius excluding the eight less regular sys- 
tems would be 213 kpc, with a median of 140 kpc. Nonethe- 
less, AMKB98 found only 6 clusters out of 60 with r core < 
100 kpc. Thus our sample spans a broader range of core radii 
(mostly due to the inclusion of less regular systems) but is 
also more peaked toward low values. This could be due to in- 
trinsic differenc es in the nature of the cl uster sample (indeed 
ISolanes. Salvador-Sole & Gonzalez-Casadol 19991 found a very low 
level of substructure in the ENACS clusters), as well as to the differ- 
ent area used, and possibly to differences in the method adopted, as 
the intrinsic coupling between r core and P in the model (we show 
for reference the r core vs. P relation in figure [121 is expected to 
produce larger core radii in correspondence of the steeper power- 
law slopes. 

On the other hand, the peak of the core radii distribution at 100 
kpc is consistent with the small values (< 100/i -1 kpc) of the core 



3 The error on /3 is ap/y/N where N is the number of clusters in the 
sample. 

4 We rescaled the AMKB98 sizes by a factor 2 since they assume a cos- 
mology with Ho = 100, go =0 



radiu s recently reported by several autho rs (e.g. lGirardi & Mezzettil 
1200 it iKatgert. Biviano & Mazurdl2004l) . Our tests (cfr. figure 
show that a systematic underestimate of our background level by 
~ 15% (or viceversa an overestimate of AMKB98) would result in 
an average r core consistent with the AMKB98 value, provided that 
the most irregular clusters discussed in § 2 are excluded from the 
sample. However the P distribution remains inconsistent at the 3a 
level. 

The fits performed on the 3x3 Mpc region are generally con- 
sistent with those performed on the whole 6x6 Mpc area, even 
though somewhat smaller, yielding a median r core — 220 kpc. 
However a KS test still rejects the hypothesis that our 3x3 Mpc 
sample and the AMKB98 one are drawn from the same distribu- 
tion at the 97% level. Note that a misplacement of the cluster center 
would flatten the central density profile resulting in larger core radii 
(see for instance iBeers & Tonrvlll986h . According to AMKB98, 
displacements smaller than 100 kpc have very small effect on the 
central profile even in the extreme case of a cusped distribution. 
Thus while such an effect could explain the few clusters with very 
large r core , it is unlikely to affect the majority of the clusters. 

Significant alignment between clu sters and their dom inant 
galax ies has often been observed (see lLambas. Groth &Peeblesl 
1988, and referenc es therein). This effect seems to be independent 
of cluster richness jFuller. West & Bridges! 1999T) and is confirmed 
up to high redshift. It involves only the cluster dominant galaxy, 
while the other bright galaxies appear to be randomly oriented, so 
that the alignment has been claimed to be p roduced by the s ame 
processes that created the dominant galaxy iKim et al]|200ll) . In 
figure fT3l we show the misalignment between the position angles 
of the clusters and those of their brightest cluster members (BCM), 
for those clusters with ellipticity large enough for the position an- 
gle to be well constrained. We note that the BCM definition itself 
depends on the observed area and that it may happen that the BCM 
lies far from the cluster centre, thus questioning the association of 
the galaxy with the cluster (if no redshift measurement is avail- 
able). Because of this reason, when the BCM was located in the 
cluster outskirts we used the position angle of the giant elliptical 
closer to the actual cluster centre, as derived from the fitting pro- 
cedure. In our case, however, this correction affects only slightly 
the alignment results, since it involves less than one third of the 
'high ellipticity sample (i.e. e > 0.25). As can be seen, the correla- 
tion existing between the two quantities is confirmed: the random 
distribution is rejected at a significance level higher than 99.9 %. 

For the other galaxies in the cluster, inspection of the right 
panels in figure [5] shows that the distribution of their position an- 
gles relative to the cluster orientation are consistent with a random 
distribution, with the possible exception of clusters A2065, A2142, 
A2178 and A2223. In fact, according to the Kolmogorov-Smirnov 
confidence levels listed in table [3] (column 9), these four clusters 
show a position angle distribution significantly different from the 
random case. 

To compare the morphological parameters with the cluster 
richness we considered two different estimates. The first (Ra) is 
calculated using a criterion similar to the Abell definition, count- 
ing how many galaxies fall in the cluster area, with magnitudes 
in the range [m.3, m.3 + 2], and then subtracting the number of 
field galaxies in the same magnitude range and in an equally large 
region. The Abell richness estimate, while being historically and 
operationally correct, is known to suffer from several drawbacks, 
mostly due to the use of a fixed apparent magnitude range (see, for 
instance lGal et all2003l) . 

For these reasons we also use a second richness estimate Rlf 
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based on the cluster luminosity function, which allows us to ex- 
ploit the entire galaxy sample down to the completeness limit. We 
used the luminosity functions (LFs) derived by Paolillo et al. (2001; 
Paolillo et al., in preparation) to compute, for each cluster, the ratio 
of the number of galaxies brighter than its absolute completeness 
limit to the number of galaxies brighter than the same limit in the 
deepest sampled cluster iGarilli. Maccagni & Andreor]l 19991) . As- 
suming that the luminosity function has a universal shape, this ratio 
provides an estimate of the richness of each cluster relative to the 
richness of the deepest sampled one. While Rlf is not exempt 
from uncertainties, in particular due to the assumption of the uni- 
versality of the luminosity function (which is challenged by recent 
data, e.g. lPiranomonte et alj 2001. Paolillo et al. in preparation), we 
note that this approach is equivalent to the use of a fixed absolute 
magnitude range, and in fact Rlf is well correlated to the richness 
estimate based on M* employed bv lGal et al 1 <2003l) . A compari- 
son between Ra and Rlf is shown in Figure fT4l 

In Figure [T5l we show the dependence of the cluster core radius 
on the cluster richness. While the core radius shows little depen- 
dence on the Abell richness, a mild correlation (correlation coeffi- 
cient of 0.62) is found using the LF richness estimate. The stronger 
correlation found with Rlf could be due to the fact that the latter 
quantity is a better estimator of the cluster richness since it takes 
into account the entire detected galaxy population, i.e. the same 
population used to measure the morphological parameters, while 
RAbeii samples only the brightest galaxy members. Furthermore, 
if the core is mainly produced by the fainter galaxies, as suggested 
by AMUS01, we expect a better correlation with Rlf than with 
RAbeii- This could also explain the lack of correlation found by 
G95, who used the Abell richness. 

The better reliability of the richness Rlf could also be sup- 
ported by the mild correlation of the cluster X-ray luminosity with 
Rlf (correlation coefficient 0.6), even with only 21 clusters for 
which X-ray luminosities are available, while no evident correla- 
tion can be seen with the Abell richness (see figure fTTl . 

Finally, we checked for correlations between the ellipticity 
and the cluster richness, as found for instance by de Theije et al. 
(1995). In figure fT6l we plot the ellipticity-richness relation for both 
Rlf and Ra', empty symbols mark the less regular clusters dis- 
cussed earlier, that are expected to show a larger ellipticity. We do 
not find any significant correlation between the ellipticity and the 
cluster richness. We note however that clusters with Ra < 40 tend 
to reach higher ellipticities, while no cluster with e > 0.6 is found 
at larger richnesses. Thus we cannot exclude that richer clusters 
tend to have, on average, smaller ellipticities. A larger sample is 
required to confirm this trend. 

No evident correlation can be seen between the cluster mor- 
phological parameters and the cluster X-ray luminosity; however, 
X-ray luminosities are available for only 21 clusters of this sample, 
and the errors on the derived morphological parameters are large, 
thus this work is not suitable to study such correlations. 



5 CONCLUSIONS 

As already stressed by many authors (see, for instance, G95), the 
derivation of cluster's morphological parameters is strongly depen- 
dent on the properties of the galaxy sample. In particular, inho- 
mogeneity of data (in terms of photometric bands, limiting magni- 
tudes, etc.), and incomplete coverage of the cluster area, may intro- 
duce systematic errors in the derived sizes, richnesses, ellipticities, 
etc. In order to minimize these problems, our sample was derived 



from a homogeneous data set extracted from DPOSS data with well 
controlled photometric errors, limiting magnitudes, and wide area 
coverage. The sample was cleaned of all objects with visually ev- 
ident signs of substructure, in an attempt to fit only compact clus- 
ters for which meaningful values of the profile parameters could 
be derived. The morphological parameters, namely the core radius, 
(3 power law index, ellipticity, position angle, and richness were 
derived, using software originally designed for X-ray data analysis, 
from the unsmoothed and essentially unbinned galaxy distributions. 

The most relevant results of this work may be summarized as 
follows: 

• One third of the clusters in this sample have core radii smaller 
than 100 kpc, which is close to the limit that our data allow us to 
resolve, and possibly consistent with cusped profiles. 

• The remaining clusters span a broad range of core radii up to 
~ 1500 kpc. While a few of these clusters are likely to be disturbed 
systems, we find several objects with regular morphology which 
seem to possess a well resolved core radius. 

• More than 80 per cent of this compact cluster sample has el- 
lipticity greater than 0.2, with an average ellipticity of ~ 0.5. 

• The comparison of the cluster position angles obtained within 
1.5 and 3 Mpc are in very good agreement, thus confirming the ab- 
sence of large substructures, and the lack of significant twisting of 
the cluster isopleths. We confirm the strong alignment between the 
cluster position angle and the major axis of the brightest ellipticals 
close to the cluster centre. 

• We find an average power-law slope of the cluster profiles 
~ 0.8). This result could be compatible with flat universe mod- 
els either with zero (standard CDM) or non-zero cosmological con- 
stant. 

• We find evidence that the core radius is correlated, albeit with 
large scatter, with the cluster richness, if the entire cluster popula- 
tion is taken into account. On the other hand, we do not find any 
significant correlation between cluster ellipticity and richness, al- 
though a mild trend cannot be excluded with our data. 

Future work will allow to confirm and extend our results, using 
the procedure developed in the present study, to estimate morpho- 
logical parameters for the larger cluster sample now available in 
the DPOSS catalogs. Furthermore, it would be desirable to apply 
the same procedure to galaxy clusters extracted from cosmologi- 
cal simulations and compare the results with real clusters, to derive 
more accurate information about the cosmological parameters and 
the process through which galaxy clusters are assembled. 
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Table 1. The compact cluster sample. The second and third columns give the cluster coordinates (from the Abell catalogue) 
at equinox J2000.0. In columns 4, 5, 6, and 7 we give, respectively, the cluster spectroscopic redshift derived from literature 
(in parentheses, the number N z of galaxies used to measure z\ see lstruble& Rood 1 19991) for discussion on lower limits), the 
Abell richness class (Abell 1958), the X-ray luminosity in the 0.1-2.4 keV band (from the BAX database - LATT), and the 
plate of the POSS archive where the cluster is imaged. 
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Table 2. Parameter estimates for the compact clusters. Column 1: cluster identification; columns 2-6: estimated parameters: 
core radius, fi power law index, ellipticity, major axis position angle and projected central density. The errors are l<r levels 
derived as discussed in the text. The tabulated So are not corrected for the different absolute completeness limits of the 
clusters. 
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Table 3. Parameter estimates for the compact clusters. Column 1: cluster identification; columns 2-3: richness estimates (Ra, 
Rlf) obtained as described in the text; column 4: total number of galaxies used for profile fitting in the 6 X 6 Mpc 2 region; 
column 5: background density as estimated in the external annulus (4 ^ r ^ 5 Mpc); column 6: background density as 
estimated with a simultaneous fit on all 8 parameters; column 7: Kolmogorov Smirnov confidence levels for the distribution 
of the cluster galaxy position angles. The errors are la levels. 
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Figure 1. Typical examples of cluster galaxy number density maps, after 
smoothing with a Gaussian filter with <r = 250 kpc (in the cluster rest- 
frame). The upper panels show two typical subclustered systems (A1035 
and A 1081) which are rejected by our selection criteria, while the lower 
panels present two compact clusters (A286 and A1661) which are included 
in our sample. Each map covers an area of 10 X 10 Mpc 2 . 
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Figure 3. Example of a unidimensional fit on binned data for one cluster of 
the sample (A286), with both a cusped profile (dotted line) and a /3-model 
(dashed line). 
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Figure 4. Effect of a 15% background misestimate on the parameters: core 
radius (top left), [3 (top right), ellipticity (bottom left) and position angle 
(bottom right). In all the panels, the diamonds mark the parameter value 
fitted with a background density overestimated by 15% of its measured 
value vs. the usual measured value, while the squares mark the parameter 
value fitted with a underestimated background vs the usual value. All the fits 
are performed in the 6 X 6Mpc 2 region. The solid line traces the bisector. 



Figure 2. Smoothed maps for A2083 (top) and A175 (bottom). Smoothing 
is gaussian with a = 250 kpc (left) and a = 150 kpc (right). The lower right 
panel shows cluster substructure on scales of order 150 kpc. The whole 
frame is 10 X lOMpc 2 wide, while the inner square is 6 X 6Mpc 2 . 
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Figure 5. Dependence of the parameters on the fitted region: core radius 
(top left), (3 (top right), ellipticity (bottom left) and position angle (bottom 
right). In all the panels, the empty/filled symbols mark the parameter value 
fitted in the 3 X 3Mpc 2 / 10 X lOMpc 2 regions respectively vs. the value 
in the 6 X 6Mpc 2 region. The dashed line traces the bisector. 
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Figure 7. Retrieved vs. true core radius for the mock cluster sample. The 
dashed line traces the bisector. 




Figure 6. Distribution of the mock cluster parameters. Upper panels: distri- 
butions of the retrieved fi values (right) and centre coordinates (left - solid 
and dotted line for xq and yo, input value was (60,60)). Lower right panel: 
distribution of retrieved ellipticity; the solid line shows the whole sample, 
dashed line denotes the high S/N clusters. Lower left panel: distribution of 
position angles, solid line is for the whole sample, dotted line is for clusters 
with retrieved ellipticity greater than 0.25. 
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Figure 8. Cluster A98. Top left panel: radial profile (the solid line marks the model while dots mark the data, see text for details). Top right panel: cluster 
isocontours in the fitted region (6x6 Mpc 2 ) plotted on the model. Bottom left and central panels: 1 and 2 a confidence regions for the couples (r CO re — /3) 
(with So free to vary) and (e — 9) respectively (r core is in units of bins (1 bin = 0.325 arcmin ~ 50 kpc), while 9 is expressed in radians). Bottom right panel: 
the misalignments of galaxy P. A. with cluster P.A. (solid line) compared with those expected for a random distribution of galaxy position angles (dashed line). 
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Figure 9. (continued). Starting from top: A171, A175, A192, A286. 
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Figure 9. (continued). Starting from top: A971, A1661, A1672, A1677. 
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Figure 9. (continued). Starting from top: A2065, A2069, A2083, A2142. 




Figure 9. (continued). Starting from top: A2177, A2178, A2223. 
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Figure 10. Distributions of the /3-model parameters for the compact cluster 
sample: ellipticity (left), /3 power-law index (center) and core radius (right). 
The solid and dashed lines show the results obtained through the fit in the 
6x6 Mpc 2 and 3x3 Mpc 2 regions, respectively. 
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Figure 11. Position angles obtained for compact clusters fitting the galaxy 
distribution within 1.5 Mpc from the cluster centre vs. position angles 
within 3 Mpc from cluster centre. The open circles refer to clusters for 
which the position angle of the major axis is poorly determined due to the 
small ellipticity (e < 0.25); the continuous line traces the bisector. 
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Figure 12. The coupling of core radius and j3 values. Different symbols 
show fit results in the three different regions. Solid dots mark the results 
from the 6x6 Mpc 2 region, while triangles and squares mark the results 
from the 3x3 and 10 X 10 Mpc 2 regions respectively. 
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Figure 13. Orientation of the cluster dominant galaxies relative to the clus- 
ter major axis (clusters for which the position angle is poorly determined 
(see text and figure lTT] are not included). 
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Figure 14. Comparison of the two richness estimates ta ( based on the 
Abell criterion) and Tl,f (based on the clusters luminosity function - see 
text for details). 
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Figure 15. Core radius versus cluster richness, for both the Abell and lumi- 
nosity function richness estimates. 
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Figure 16. Ellipticity vs. cluster richness, for both the Abell and luminosity 
function richness estimates. Empty symbols mark the less regular clusters 
(see text for details). 
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Figure 17. Cluster X-ray luminosity vs. richness, for both the Abell and 
luminosity function richness estimates. Only 21 clusters with X-ray lumi- 
nosities available are shown . 



